Supercritical bifurcation of a hula hoop 
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The motion of a hoop hung on a spinning wire provides an illustrative and pedagogical example of a supercrit- 
ical bifurcation. Above a certain angular velocity threshold Q, c , the hoop rises, making an angle 6 ~ (SI — f2 c ) 2 
with the vertical. The equation of motion is derived in the limit of a long massless wire, and the calculated steady 
states are compared to experimental measurements. This simple experiment is suitable for classroom demon- 
stration, and provides an interesting alternative to the classical experiment of the bead sliding on a rotation 
hoop. 



I. INTRODUCTION 

The rotation of rigid bodies often displays interesting insta- 
bility problems. Bodies with three different inertia moments 
are well known to have unstable rotation about the interme- 
diate axistl, as commonly observed from acrobatic jumps or 
dives. Throwing a tennis racket provides an easy illustration 
of this instability of the free rotation. Such purely inertial 
instability has no threshold, i.e. it can be observed for very 
low rotation rates. In contrast, for axisymmetric bodies such 
as hoops, disks, rockets, eggs etc., the free rotations remain 
stable, but different mechanisms may also lead to instability 
when external forces are present. Under some circumstances, 
spectacular and unexpected instabilities may originate from 
frictional forces. This is the case for the tippe-top, a popular 
toy that flips over and rotates on its stem, or for the hard- 
boiled egg problem, which has recently received a nice anal- 
ysisa. More classically, the competition between gravity and 
centrifugal forces may lead to an instability with a finite ro- 
tation rate thresholdO, as illustrated by the simple experiment 
described in this paper. 

A hoop is hung on a long wire, whose upper end is spun. At 
low rotation rate, the hoop is vertical and simply spins about 
its diameter. Increasing the rotation rate, the hoop progres- 
sively rises and becomes horizontal, spinning about its sym- 
metry axis. This situation may appear paradoxical, since the 
horizontal position maximizes both kinetic and potential en- 
ergy. It is similar to the conical pendulum problerro, com- 
monly illustrated in the classical demonstration.axperiment of 
the bead sliding along a vertically rotating hoopcBtl. The pop- 
ular 'hula hoop' game, where the wire rotation is replaced by 
the hips oscillations of the player, is a common illustration of 
this phenomenon. Once the hoop spins horizontally, its rota- 
tion is maintained by a parametric oscillation mechanism!]. 

In this paper, the equation of motion is derived from the La- 
grange's equation. An alternate derivation, from the angular 
momentum equation, is also presented as a good illustration 
of dynamics of a rigid body. The hoop is shown to rise fol- 
lowing a supercritical bifurcation for the wire rotation rate 57 
above a critical value Jl c = {2g/ R) 1 / 2 , where g is the gravita- 
tional acceleration and R the hoop radius. The stable solution 
6 ~ (f2 — 51c) 1 / 2 coincides with that of the conical pendulum 
or the bead-on-a-hoop. The nonlinear oscillations are briefly 
described by means of phase portraits, and compared to that of 
the bead-on-a-hoop. Simple considerations from wire torsion 




FIG. 1: A hoop of center of mass G, attached to a wire OP, spins 
about the vertical axis OG, making a pitch angle 9 with the vertical. 



and air friction allow us to estimate the startup and damping 
timescales. Finally some experimental measurements are re- 
ported, and are shown to compare well with the exact solution. 

In this experiment, the angular rotation threshold being of 
order of a few rad/s, the wire can be simply spun with the 
fingers. Since it only requires a wire and a hoop, this exper- 
iment can be conveniently used as a classroom illustration of 
spontaneous symmetry breaking and bifurcation. Although 
the detail of the calculation is somewhat more subtle than that 
of the classical bead-on-a-hoop problem, since it deals with 
rigid body dynamics, the physics is basically the same and 
does not require the much heavier apparatus of the bead-on-a- 
hoop demonstration. 
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II. THEORY 
A. Equation of motion 

We consider a thin uniform hoop, of radius R and mass m, 
fixed at a point P of its periphery by a massless wire of length 
L ^> R (see figure |]). The wire is spun from its other end 
O at a constant angular velocity £7. The center of mass G is 
assumed to remain at the vertical of O. Let us first consider a 
torsionless wire, so that the hoop angular velocity <j> about the 
vertical axis e 2 instantaneously follows the imposed angular 
velocity f2. 

Let us consider the rotating frame of reference of the hoop 
(ei , e2 , e3 ) . We note 9 the pitch angle between the hoop plane 
and the vertical axis e z . We consider for simplicity a zero 
roll angle about e! (i.e. we assume that the point P remains 
the highest point of the hoop), so that 9 is the only degree of 
freedom for this problem. The absolute angular velocity w of 
the hoop has two contributions. The first one, of magnitude 57, 
is imposed from O through the torsionless wire, and is about 
e z = cos6>ei + singes (the points O and G are at rest, so 
that OG is the instantaneous axis of rotation when 9 is kept 
constant). The second one comes from the pitch variation 9, 
and is about the axis e2. We therefore obtain 



(1) 



The potential energy is V = mgzQ, where zq ~ R(l — 
cos 9) is the center of mass elevation to first order in R/L. The 
kinetic energy T has two contributions: one from the rotation 
■ I ■ lj, where / is the inertia matrix of the hoop relative 
to G, and one from the vertical translation of the center of 
mass, \mzQ ~ ^mR 2 9 sin 2 9. In the reference frame of the 

hoop, the inertia matrix / has diagonal components h = I2 = 
\mR 2 and I3 = mR 2 . The Lagrangian function C = T — V 
then finally writes 




1 



C = -m,R 2 
4 



51 2 (l + sin 2 9) + 9 2 {l + 2 sin 2 



-mgR cos 9. 
(2) 

Writing the Lagrange's equation for the coordinate 9, we end 
up with the differential equation of motion 

0(1 + 2 sin 2 9) - (Q 2 - 29 2 ) sin 9 cos 9 + Q 2 C sin 9 = 0. (3) 



where Q c — \j2gjR. 

This equation of motion can also be obtained from the an- 
gular momentum equation and Newton's law. In the rotating 
frame of reference, the angular momentum equation reads 



^ + " XL = r ' 



(4) 



where L = Iu> is the angular momentum and F = GP x N is 
the torque of the wire tension relative to the center of mass (the 
gravitational torque vanishes). The wire tension is obtained 
from Newton's law, 



mz G e z 



N 



TOg. 



(5) 



To first order in R/L, N is vertical, so its torque is simply 
r ~ —RNs'm9e2- Using zg — R(l — cos#), we obtain 



-R (mg + mR(9 2 cos 9 + 9 sin 9)j sin 9 e 2 . 



(6) 



Replacing into (|j) and projecting on e2, we finally recover the 
equation of motion (||). 



B. Steady states and small oscillations 

We are interested in the steady states and the natural fre- 
quency for small oscillations about them. To first order in 9, 
equation O) reduces to 



9 - ST sinf 



nj sin 9 = 0. 



(7) 



It is worth noting that this linearized equation is exactly the 
same as the one from the bead-in-a-hoop problemD. The dif- 
ference between equations (||) and (0) originates from the ki- 
netic energy of the center of mass translation, which is not 
present in the bead-on-a-hoop. 

In addition to the trivial solution 9 cq = 0, equation (Q) has 
a non trivial solution for > fL, 



±COS [- 



(8) 



Following the usual terminology, the pitch angle 9 eq is the 
order parameter, and we introduce the reduced angular veloc- 
ity e = il/il c — 1 as the control parameter. The lineariza- 
tion of (|8|) finally leads to the classical form of a supercritical 
pitchfork bifurcation for e«l: 



(9) 



Stability and natural frequency for small oscillations are ob- 
tained by introducing in equation (0) a small perturbation in 
the form 



9(t) = 



= (cr + iw)t 



(10) 



where 9 cq stands for the trivial or non-trivial solution. The 
trivial solution is stable for f2 < f2 c (a — 0), with a natural 
frequency lj — (fl 2 — 57 2 ) 1 / 2 , and unstable for Q > Q c , with 
a growth rate given by a = (Q 2 — 5! 2 ) 1 / 2 (and uj = 0). The 
non-trivial solution (||) is found to be always stable, with the 
same natural frequency. In terms of e, one can see that the 
period of the oscillations, 



T 



1-1/2 



(11) 



diverges as one approaches the transition from both side. This 
critical slewing down is a usual signature of supercritical 
bifurcatioro 



C. Nonlinear oscillations 

When the oscillation amplitudes about the steady states are 
not small, the nonlinear terms in (Jfy become important and 
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FIG. 2: Orbits in the phase space [0, 9/£l c ) for the hula hoop (solid 
lines) compared to the bead-on-a-hoop system (dashed lines), (a): 

n = o.8ti c . (by. n = i.2n c . 



interested in the closed orbits in the domain \9\ < tt/2: orbits 
that cross 9 = ±7r/2 are not consistent with the assumption 
that the point P remains the highest point of the hoop. 

Around the stable fixed points, both the two systems show 
nearly elliptic orbits, as expected from small harmonic oscil- 
lations. Larger oscillations of the hula hoop show orbits with 
sharper corners around 9 ~ 0, associated to larger angular 
velocities \8\. This discrepancy originates from the vertical 
translation of the center of mass, which is responsible for an 
additional pitch angle acceleration when the hoop is nearly 
vertical. 



D. Timescales 

Two timescales are relevant for a practical experiment: the 
startup timescale of rotation r s and the damping timescale of 
oscillations ry. The startup timescale may be obtained consid- 
ering the wire torsion. We start from an initially vertical hoop 
(6 = 0), and let now its rotation angle <fi about e z be free. 
The wire communicates a torque k{4> — fit) = nScf) to the 
hoop, where k is the torsion constant and 5<fi the wire torsion. 
The angular momentum equation about e z applied to the hoop 
then reads 

^mR 2 54> + K(50 = O, (13) 

leading to a startup timescale t s ~ R(m/2n) 1 / 2 (which cor- 
responds also to the period of the wire torsion oscillations if 
no damping were present). 

The damping timescale is relevant both for the <p oscilla- 
tions due to the wire torsion, and for the 9 oscillations around 
the equilibrium states. It can be obtained considering the dis- 
sipation with the surrounding air. With typical velocity of or- 
der m/s and hoop thickness w of a few millimeters, we can 
make use of a turbulent estimate for the drag force. Neglect- 
ing 9 compared to <f>, a unit length dl of the hoop experiences 
a drag force df ~ pR 2 <p 2 w dl, where p is the air density. Inte- 
grating over the hoop perimeter leads to the frictional torque 



may affect the dynamics of the system. The orbits in the phase 
space (9, 9) then provide a useful tool to characterize the non- 
linear dynamics of the hula hoop, and to compare it with the 
one of the bead-on-a-hoop (Q). The equation of motion can be 
integrated using the Painleve invariant, 



C = 9^ 
89 



C, 



(12) 



which takes a constant value along an orbit (this invariant cor- 
responds to the total energy E = T + V in the case of a 
conservative system). Figure || shows isolines of C for two 
forcing frequencies, below (fi/fi c = 0.8, fig. a) and above 
(fi/fic = 1.2, fig. b) the transition, together with the orbits of 
the bead-on-a-hoop problem (dashed lines). One can clearly 
see from these phase portraits the trivial solution 9 cq = for 
fl/fl c — 0.8 and the two non trivial solutions (^]) separated by 
a saddle point at 9 = for 51/57 c = 1.2. Note that we are only 



The angular momentum equation about e z reads 



pow 
p R 



where po is the hoop density. Taking 
the damping timescale finally writes 



£1 and 



Td 



p R 



(14) 
(15) 

(16) 



For values of practical interest, is of order of a few tens 
of the rotation period O -1 . Note here that in the case of the 
real hula hoop game, the rolling friction on the player hips is 
fortunately dominant, leading to a damping timescale of order 
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III. EXPERIMENTAL RESULTS 



(a) 




(c) I (d) 



FIG. 3: Pictures of the hoop, (a) The hoop at rest, (b) and (c) 9 ~ 
37°. (d) 9 ~ 78°. 

An experiment has been carried out using a wood hoop 
(po ~ 0.67 g.cm~ 3 ) of mass m = 15.5 g, radius R = 86 mm 
and section 4 x 10 mm 2 . The expected frequency threshold is 
then f2 c /27r ~ 2.40 Hz (classical hula hoops have H c /2tt ~ 
1 Hz). A simple cotton thread, of length L = 0.68 m and mass 
less than 3% of the hoop mass, was fixed to a constant current 
motor. The startup timescale for this wire can be estimated 
from the free oscillation period, r s ~ 5 s (corresponding to a 
torsion constant k ~ 2 x 10~ 6 N.m). This wire is far from 
being torsionless, and during the early stage of the rotation, 
torsional energy is stored into the wire. The wire progres- 
sively communicates rotation to the hoop, and the fluctations 
(due to a varying imposed rotation rate at its upper end, or to 
pitch angle variation at the other end) are smoothed down on 
a timescale r^. 

Seen from the side, the hoop rotation appears as Lissajous 
ellipses, whose principal axis makes an angle 9 with the ver- 
tical axis. Pictures acquired from a simple CCD camera (see 
figure |J) allows us to easily measure the pitch angle 9 within 
1° from these ellipses. The time aperture of the camera may 
blur the picture, but facilitate in some case the angle measure- 
ment (see for instance picture c). 

The measured pitch angles 9, shown in figure Q are in 
excellent agreement with the exact solution (|8|), except for 
low rotation rate, where non zero angles are measured below 
the expected transition. The experimental frequency thresh- 
old, obtained by extrapolating the curve down to 9 ~ 0, is 
fl c /2TT ~ (2.30 ± 0.03) Hz, which agrees within 5% with the 
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FIG. 4: Bifurcation diagram of the pitch angle 9 as a function of the 
normalized angular velocity Q/Q c . The solid line is the calculated 
solution (p|), with Q c adjusted to fit the experimental data. 



theoretical value. The discrepancy at low rotation rate is typi- 
cal of an imperfect bifurcation, where small asymmetry in the 
apparatus (such as the position of the knot) slightly anticipates 
the destabilization of the basic state. 

An interesting observation is that, for fl > 2.6il c (see the 
dashed line in figure Q), the bifurcated state is no stable any 
more: a secondary instability appears, in the form of a slow 
precession of the center of mass, with a period of about ten ro- 
tation periods. This new behavior is probably an effect of the 
non zero mass of the wire, and can clearly not be described in 
our calculation, where the hoop center of mass is constrained 
to remain aligned with the vertical axis. 

Transient phenomena may also be investigated, by means 
of usual video processing. An illustration is given in figure |[ 
showing a spatio-temporal diagram obtained by collecting the 
light intensity recorded on a horizontal line passing through 
the wire. In this example, the imposed frequency Q/2w has 
been suddenly increased from to 5.5 Hz ~ 2.4 £l c /2ir (first 
arrow). One can see that, after a transient time of around 
7 s (second arrow) during which the wire rotation propagates 
down to the hoop, the amplitude increases and saturates to a 
finite value. 

From this diagram, the instantaneous pitch angle as well 
as the instantaneous oscillation frequency (f>/2ir may be ex- 
tracted, as shown in figure |[ The frequency is obtained from 
averaging over 6 successive oscillations. Due to the progres- 
sive torsion of the wire, this frequency slowly approaches its 
imposed value 5.5 Hz. As a consequence, the increase of 
the pitch angle towards its stationary value 9 cq ~ 1.42 rad 
is rather slow. 

It is interesting to note that, when plotting the instantaneous 
pitch angle as a function of the instantaneous frequency (see 
figure Q), the bifurcation diagram of figure ^ is recovered to 
a high degree of accuracy. This suggests that the pitch an- 
gle follows 'adiabatically' the instantaneous frequency, so that 
the apparent growth rate is essentially controlled by the wire 
torsion more than the intrinsic dynamics of the instability (at 
least far from the transition). 
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FIG. 5: Spatio-temporal diagram of the light intensity recorded on 
an horizontal line, showing the wire oscillation, when the frequency 
Q/2n is increased from to 5.5 Hz. 




FIG. 6: Instantaneous pitch angle 9 (left) and rotation frequency 
<j)/2n (right) as measured from figure pi The origin time to corre- 
sponds to the second arrow in that figure. 



IV. DISCUSSION 

This simple experiment provides an interesting alternative 
to the bead-on-a-hoopo3tl experiment as a mechanical analog 
to bifurcation and second-order phase transition in physics. 
The poor attention paid to the choice of the material and the 
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0.5 1 1.5 2 2.5 3 
d>/Q 

T c 

FIG. 7: Instantaneous pitch angle 8 as a function of the instantaneous 
rotation frequency <j>/Q c (same data as in figure ^). The solid line is 
the calculated solution (^|), as in figure ^. The cross indicates the 
stationary solution # cq ~ 1.42 rad for Q — 2.4 Q c . 



experimental conditions illustrates the robustness of the phe- 
nomenon, and makes this experiment easy to work out by un- 
dergraduate students. Further experiments can be performed, 
e.g. studying the transient phenomena resulting from a small 
perturbation. As a suggestion, restrain the point P on the ver- 
tical axis when £1 > f2 c by means of a small hook around the 
wire. Releasing the hook allows us to measure the growth rate 
and to characterize its divergence as e — * 0. This method can 
be hardly achieved with a bead-on-a-hoop apparatus due to its 
inherent difficult access. 

Another motion may compete with the bifurcation de- 
scribed in this paper: the hoop can rotate as a whole about 
the vertical axis, its center of mass remaining aligned with 
the wire. This is the usual motion for the conical pendulum, 
which also leads to a supercritical bifurcation with a threshold 
simply given by the natural frequency luq ~ \fg] (L + R). 
The conical pendulum motion overcomes the hula hoop mo- 
tion for luq of order of il c , i.e. for a wire length of order of 
the hoop radius. Moving away the center of mass from the 
vertical axis may allow us to observe the competition between 
the two regimes, even for a longer wire. 

Similar bifurcations as the result of a competition between 
centrifugal force and gravity are present in a number of sit- 
uations. Spinning plates provides an interesting illustration: 
as for the real hula hoop, the motion here is forced by the 
precession of the rigid rod rather than its rotation. The im- 
posed precession frequency has to overcome the natural fre- 
quency for the plate to stand up, and then the horizontal state 
is maintained by a parametric oscillation mechanism!! Lasso 
roping is another example, where the rigid hoop is replaced 
by a deformable loop. Here again the motion is maintained by 
the precession of the spoke rather than rotation. As a conse- 
quence, the knot on the loop makes the spoke to rotate as well, 
so that the cord has to be continuously untwisted at its other 
end. 

It is worth pointing out that this experiment can be car- 
ried out with any rigid body, not necessarily axisymmetric. 
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In this case, the angle of the equilibrium state (||) remains 
unchanged, and the angular velocity threshold just becomes 
fl c = y/mga/ (I3 — where a = | GP | is the distance be- 
tween the center of mass and the knot. Note however that, 
for bodies with three different inertia moments, in addition to 
the bifurcation described here, inertial instabilities may also 
occur. Such instability involves the roll angle about ei as an 
additional degree of freedom, and is not described by our cal- 



culation. 
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